Coherent backscattering effect on wave dynamics in a random medium 
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A dynamical effect of coherent backscattering is predicted theoretically and supported by computer 
simulations: The distribution of single-mode delay times of waves reflected by a disordered waveguide 
depends on whether the incident and detected modes are the same or not. The change amounts to 
a rescaling of the distribution by a factor close to y/2. This effect appears only if the length of the 
waveguide exceeds the localization length; there is no effect of coherent backscattering on the delay 
times in the diffusive regime. 



PACS numbers: 42.25.Dd, 42.25.Hz, 72.15.Rn 

Coherent backscattering refers to the systematic con- 
structive interference of waves reflected from a medium 
with randomly located scatterers. The constructive inter- 
ference occurs in a narrow cone around the angle of inci- 
dence, and is a fundamental consequence of time-reversal 
symmetry jlj]. The resulting peak in the angular depen- 
dence of the reflected intensity is a generic wave effect: 
It has been observed using light waves || and acoustic 
waves H , for classical and quantum scatterers Q , in pas- 
sive and active media Q . 

These studies mainly addressed static properties. Dy- 
namic aspects of wave propagation in random media are 
now entering the focus of attention |^-[|, and the work 
on acoustic waves || has started to study the connec- 
tion with the coherent backscattering effect. The key 
observable in the dynamic experiments Q is the deriva- 
tive (j)' = d(f>/duj of the phase <fi of the wave amplitude 
with respect to the frequency lu. The quantity <j)' has the 
dimension of a time and is interpreted as a delay time. 
Van Tiggelen et al. M have developed a statistical theory 
for the distribution of the delay time <f>' and the intensity 
J in a waveguide geometry (where angles of incidence are 
discretized as modes). Although the theory was worked 
out mainly for the case of transmission, the implications 
for reflection are that the distribution P (</>') does not de- 
pend on whether the detected mode n is the same as the 
incident mode m or not. This is in contrast with P(I), 
which is rescaled by a factor of 1 /2 when n becomes equal 
to to — so that the mean / becomes twice as large. Hence 
it appears that no coherent backscattering effect exists 
for P{<j>'). 

What we will demonstrate here is that this is true only 
if wave localization may be disregarded. Previous studies 
HQ dealt with the diffusive regime of waveguide lengths 
L below the localization length £. Here we consider the 
localized regime L > £ (assuming that also the absorp- 
tion length £ a > £). The distribution of reflected inten- 
sity is insensitive to the presence or absence of localiza- 
tion, being given in both regimes by Rayleigh's law: 
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(for unit incident intensity). In contrast, we find that 
the delay-time distribution changes markedly as one en- 



ters the localized regime, decaying more slowly for large 
\<fr'\. Moreover, a coherent backscattering effect appears: 
For L > £ the peak of P{4>') is higher for n — to than 
for n ^ to by a factor which is close to \/2. We present 
a complete analytical theory, compare it with numeri- 
cal simulations, and offer a qualitative argument for this 
unexpected dynamical effect of coherent backscattering. 

Let us begin with a more precise formulation of the 
problem. We consider a disordered medium (mean free 
path I) in a waveguide geometry (length L, with N 3> 1 
propagating modes at frequency u>) and study the cor- 
relator p of the reflected wave amplitudes at two nearby 
frequencies ui ± \5lo, 



p = r nm (uj + Uuj)r* nm (uj - Uuj) 
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The indices n and to specify the detected and incident 
mode, respectively. (We assume single-mode excitation 



and detection.) The amplitudes 



form the N x N 



reflection matrix r. In the localized regime (localization 
length £ ~ Nl smaller than both L and the absorption 
length £ a ), the matrix r is approximately unitary because 
transmission is negligibly small. We assume time-reversal 
symmetry (no magneto-optical effects), so that r is also 
symmetric. Following Genack etal. [g|J^] we define the 
single-mode (or single-channel) delay time </>' as 
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where I = \r nm (uu)\ 2 is the intensity of the reflected wave 
in the detected mode for unit incident intensity. We seek 
the joint distribution function P(I, (j)') m & n ensemble of 
different realizations of disorder. 

The single-mode delay time <fi' is a linear combination 
of the Wigner-Smith jlQ] delay times Tj (i = 1, 2, . . . , N), 
which are the eigenvalues of the matrix 
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(The matrix of eigenvectors U is unitary for a unitary 
reflection matrix.) For small Slo we can expand 

r(w ± Uuj) = U T U ± hi Slu U t diag (t u ...,t n )U, (5) 
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hence the relations 
At 



y = Re^i, I=\A \ 2 , A k = Y j r?u i v i . (6) 



We have abbreviated Ui — Ui m , V{ — Ui n . 

The distribution of the Wigner-Smith delay times for 
this problem was determined recently Q. In terms of 
the rates fii ~ 1 /r, it has the form of the Laguerre en- 
semble of random-matrix theory, 
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where Q(x) = 1 for x > and for x < 0. The pa- 
rameter j — al/c (with wave velocity c) equals the scat- 
tering rate times a numerical coefficient (a — 7r 2 /4, 8/3 
for two, three-dimensional scattering). Eq. (|7|) extends 
the single-mode (N = 1) result of Refs. [p~2|— |l4[| to any 
N. The matrix U is uniformly distributed in the uni- 
tary group. We consider first the typical case n ^ m of 
different incident and detected modes. (The special case 
n = m is addressed later.) For n ^ m the vectors u and 
v become uncorrelated in the large- N limit, and their el- 
ements become independent Gaussian random numbers 
with vanishing mean and variance {\uf |) = |) = AT . 

It is convenient to work momentarily with the weighted 
delay time W — <j>'I and to recover P(I, <fi') from 
P{I, W) at the end. The characteristic function \{Pi q) = 
(e-ipi-w~W) is the Fourier transform of P(I,W). The 
average (• • •) is over the vectors u and v and over the 
set of eigenvalues {n}. The average over one of the vec- 
tors, say v, is easily carried out, because it is a Gaussian 
integration. The result is a determinant, 
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The Hermitian matrix H is a sum of dyadic products 
of the vectors u and u, with Ui = UiTi, and hence has 
only two non- vanishing eigenvalues A+ and A_. Some 
straightforward linear algebra gives 



\± = \{qB 1 +p± ^2pqB 1 + q 2 B 2 



where we have defined the spectral moments 
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The resulting determinant is det(l + iJ/A") -1 
X+fN)- 1 ^ + \_/N)-\ hence 
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An inverse Fourier transform, followed by a change of 
variables from I,W to I, (/)', gives 



P{I,(j)') = e{I)(N 3 I/TT)^ 2 e- NI 
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The average is over the spectral moments B\ and B 2l 
which depend on the Mi's and Tj's via Eq. (11). 

This result in the localized regime is to be compared 
with the result of diffusion theory MM , 



? diff (/,^') = e(/)(Jv 3 /A) 
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The constants are given by Q ~ L/l and <j)' ~ L/c up to 
numerical coefficients of order unity jl5| . Comparison of 
Eqs. (|l3| ) and ([TJ) shows that the two distributions would 
be identical if statistical fluctuations in the spectral mo- 
ments £>i, B 2 could be ignored. However, as we shall 
see shortly, the distribution P(Bi,B 2 ) is very broad, so 
that fluctuations can not be ignored. The large fluctu- 
ations are a consequence of the high density of anoma- 
lously large Wigner-Smith delay times r t in the Laguerre 
ensemble (^), and are related to the penetration of the 
wave deep into the localized regions. The large Tj's are 
eliminated in the diffusive regime L < £. Then Bi and 
B 2 can be replaced by their ensemble averages, and the 
Gaussian theory (|J^] is recovered. (The same applies if 
the absorption length £ a < £.) 

To determine how the statistical fluctuations in the 
spectral moments alter P(I, <$>'), we need the joint dis- 
tribution P{B-y, B 2 ). This can be calculated by applying 
the random-matrix technique of Refs. [ p~6|JT^ ] to the La- 
guerre ensemble. The result is 
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where Ei (x) is the exponential-integral function. The 
most probable values are Bi ~ 7 AT, B 2 ~ 7 2 Af 3 , while 
the mean values (-Bi), (B 2 ) diverge — demonstrating the 
presence of large fluctuations. The distribution P(I, </>') 
follows from Eq. ([[3]) by integrating over B\ and B 2 with 
weight given by Eq. ([l5|). This is an exact result in the 
large- N limit. 

For the discussion we concentrate on the distribution 
P(4>') = J °° dJ P(I, 4>') of the single-mode delay time by 
itself, which takes the form 



dSi dB 2 







P{B U B 2 )(B 2 -B 2 ) 
2{B 2 + (j> 12 - 2Bi0') 3/2 



(16) 



2 



X 

K 



X 



-o- 

K 




X 



0.006 



0.004 



0.002 



0.000 



0.002 



0.001 



0.000 



-200 -100 100 200 




-600 -400 -200 200 400 600 
f /y 

FIG. 1. Distribution of the single-mode delay time <j>' in 
the diffusive regime (top panel), intermediate regime (mid- 
dle panel), and localized regime (lower panel). The results of 
numerical simulations (data points) are compared to the pre- 
diction ( |r?| ) of diffusion theory (dashed curve) and the 
prediction (jl6|) for the localized regime (solid curve). These 
are results for different incident and detected modes n 7^ m. 

We compare this distribution in the localized regime with 
the result of diffusion theory 



Pdiff(f) = (Q/20')[Q + (W-i) 
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In the localized regime the value <^ poak — 7^ at the 
center of the peak of P(<fi') is much smaller than the 
width of the peak A(f>' ~ 7 iV 3 / 2 ~ <t>' peak (£/l) 1/2 - This 
holds also in the diffusive regime, where 4>' pcak — <f>' and 
A<f>' ~ 0p eak (i/O 1/2 - However, the mean (0') = {B x ) 
diverges for P, but is finite (equal to 4>') for Pdiff- I 11 the 
tails P decays oc |</>'|~ 2 , while Pdiff oc |</>'|~ 3 - 

The transition from the diffusive to the localized 
regime with increasing L is illustrated in Fig. [j]. The 



data points are obtained from the numerical simulation of 
scattering of a scalar wave by a two-dimensional random 
medium. The reflection matrices r(ui ± \5tu) are com- 
puted by applying the method of recursive Green func- 
tions |L8| to the Helmholtz equation on a square lattice 
(lattice constant a). The width W = 100a and the fre- 
quency u> = 1.4 c/a are chosen such that there are N = 50 
propagating modes. The mean free path I — 14.0 a is 
found from the formula T — (1 + s) -1 for the transmis- 
sion probability in the diffusive regime, where s = 2L/irl 
for two-dimensional scattering. The corresponding lo- 
calization length £ = NL/s = 1100 a. The parameter 
7 = 46.3 a/c is found from </>' in the diffusive regime p9| . 
The relationship between the parameters 7, 4>' , and Q 
appearing in P and Pdiff is given by fl5| 
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In Fig. [j], the same set of parameters is used for all 
lengths to plot the distributions P (solid curve) and Pdiff 
(dashed). The numerical data agrees very well with the 
analytical predictions in their respective regimes of va- 
lidity. 

We now turn to the case n = m of equal-mode excita- 
tion and detection. The vectors u and v in Eq. (^) are 
then identical, and we can write 
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The joint distribution function of the complex num- 
bers Co and C\ can be calculated in the same way as 
P(Bx,B 2 ). We find 



P(C ,Ci) ocexp(-iV|C , o|72) / dxx 2 e~ x 

Jo 

I 1 + - -^7 ReC ci] 2 . (20) 
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The maximal value P( 



. Apeak) = yJV* N3 l 2 for n = m is 
larger than the maximum of P{4>) for n 7^ to by a factor 
y/2x 1371^ = 1-35 in the large- N limit. This is in contrast 
to the diffusive regime, where there is no difference in the 
distributions of single-mode delay times for n = m and 
n 7^ to. Our analytical expectations are again in excellent 
agreement with the numerical simulations, presented in 
Fig.| 

In order to explain the coherent backscattering en- 
hancement of the peak of P{4>') in qualitative terms, 
we compare Eq. ( |l9| ) for n — to with the correspond- 
ing relation @ for n 7^ to. The quantities Aq and 
Ai, as well as the quantities Co and C\, become mu- 
tually independent in the large- TV limit. [The cross-term 
( 7 7V)- 1 ReC C 1 * in Eq. @ is of order N' 1 ' 2 because 
C ~ N- 1 ' 2 and C x ~ lN.} The main contribution to 
the enhancement of the peak height, namely the factor 
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FIG. 2. Same as Fig. |l|, but now comparing the case 
n/mof different incident and detected modes (solid circles) 
with the equal-mode case n — m (open circles). The curve 
for n = min the lower panel is calculated from Eqs. (jl9|) and 

©■ 



of y/2, has the same origin as the factor-of-two enhance- 
ment of the mean intensity I. More precisely, the rela- 
tion P(A ) = V2P{V2C ) leads to a rescaling of P(J) 
for n — ui by a factor of 1/2 [see Eq. ([!])] and to a rescal- 
ing of P(4>') by a factor of The remaining factor of 
.tl 9 , 6 = 0.95 comes from the difference in the distribu- 
tions P(Ai) and P(Ci). These distributions turn out to 
be very similar, hence the factor is close to unity. The 
asymptotic independence of Aq and A\ (as well as of Co 
and C\) is another consequence of the strong fluctuations 
originating from the high density of anomalously large 
Wigner-Smith delay times r, . In the diffusive regime the 
corresponding quantities are strongly correlated, and the 
coherent backscattering enhancement of the intensity af- 
fects both in the same way. Because only their ratio 
features in 0', this effect cancels and no difference is ob- 
served in Pdiff (</>') for n = m and n =/= m. 

In conclusion, we have discovered a dynamical effect of 
coherent backscattering, that requires localization for its 
existence. Computer simulations confirm our prediction, 
which now awaits experimental observation. 
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